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ABSTRACT 


The study: (1) Investigates the effects of a uniform temperature change 
on the stresses and deformations of composite tubes; (2) Determines the 
accuracy of an approximate solution based on the principle of complementary 
virtual work. Interest centers on tube response away from the ends and so a 
planar elasticity approach is used. For the approximate solution a piecewise 
linear variation of stresses with the radial coordinate Is assumed. The 
results from the approximate solution are compared with the elasticity 
solution. The stress predictions agree well, particularly peak interlaminar 
stresses. Suprlsingly, the axial deformations also agree well. This despite 
the fact that the deformations predicted by the approximate solution do not 
satisfy the interface displacement continuity conditions required by the 
elasticity solution. The study shows that the axial thermal expansion 
coefficient of tubes with a specific number of axial and circumferential 
layers depends on the stacking sequence. This Is In contrast to classical 
lamination theory which predicts the expansion to be independent of the 
stacking arrangement. As expected, the sign and magnitude of the peak 
Interlaminar stresses depends on stacking sequence. For tubes with a specific 
number of axial and circumferential layers, thermally-induced interlaminar 
stresses can be controlled by altering stacking arrangement. 


INTRODUCTION 


Tubes are efficient structural elements for reacting bending, 
axial, and torsional loads. In addition, when using fiber-reinforced 
materials, tubes can be efficiently made using filament winding or one 
of several other automated techniques. Stiffness, strength, and thermal 
expansion properties of the tube can be controlled by varying the number 
of layers of fiber-reinforced material, and by varying the angle of the 
fibers relative to the axial direction. However, for fiber-reinforced 
materials like graphite-epoxy, tubes are made by curing the epoxy at an 
elevated temperature. Typical epoxies cure at 120- 175°C (250-350°F) 
while epoxies made for high temperature operation cure at 290-340°C 
(550-650°F). Operating a graphite-epoxy tube below the cure temperature 
of the material can result in high residual stresses. If the stresses 
are too high, material failure can occur, resulting In a loss of 
performance of the tube. 

Graphite-epoxy tubes are currently being considered as a major 
building block for large orbiting space structures. Ambient 
temperatures as low as -156°C (-250°F) can occur in space. To use tubes 
effectively, understanding their response In these thermal environments 
is quite important. This paper reports on the calculation of the 
response of layered tubes to a spatially uniform temperature change. 

The change represents the difference between the tube fabrication 
temperature and the temperature of its operating environment. The 
responses of interest in this paper are the thermally induced stresses 
and the axial thermal expansion characteristics. The latter response is 
important for dimensional stability considerations, both during the 
assembly of large structures, and during operation. Unless dimensional 
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changes are controlled, orbiting Into and out of the earth's shadow can 
cause problems. Of particular concern Is the Influence of layering 
arrangement on stresses and axial expansion. 

The tubes considered are of a constant Inner and outer diameters. 
The tubes are made of multiple layers, each layer having properties 
which represent unidirectional graphite-epoxy. Each layer will be 
assumed to be homogeneous and either orthotropic or transversely 
Isotropic in a coordinate system aligned with the axis of the tube. In 
terms of fiber-reinforced material, the study is restricted to what are 
referred to as cross-ply tubes. A layer is said to be a 0° layer if the 
stiff, or fiber, direction is axial. A layer is said to be a 90° layer 
if the fiber direction Is circumferential . Stacking sequence notation 
common to the analysis of layered composite materials will be used. A 
tube designated as a (90/0/90/0) tube is a cross-ply tube and has four 
layers. The inner layer and the second layer from the outside have 
their fibers aligned circumferentially (90°). The outside layer and the 
second layer from the inside have their fibers aligned axially (0°). A 
(0/90/0/90) tube would have the fibers In the Inner layer oriented 
axially, etc. Interest centers on the portions of the tube away from 
the ends. In practice, the end fittings used to connect the tubes 
together will dictate to a large degree the stress state near the 
ends. Since these regions on the end of the tube constitute such a 
small percentage of the total tube length, the overall axial thermal 
expansion characteristics, for example, could well be determined by the 
large percentage of the length of the tube away from the ends. 

Therefore, a generalized plane deformation elasticity solution is used 
to determine the tube response. The solution assumes that the stresses 
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are not a function of the axial coordinate of the tube. An elasticity 
solution, as opposed to a shell-like approach, Is used because some 
tubes being considered for use have mean radius to wall thickness ratios 
of less than 10. For such a small ratio, through-the-thlckness effects 
have to be accounted for. To account for these effects, a higher-order 
shell theory could be used. However, an elasticity approach Is felt to 
be more direct. 

With a spatially uniform temperature, the problem is assumed to be 
axi symmetric. Thus, with no dependence of the stresses on the axial 
coordinate, and using the axlsymmetry assumption, the equations 
governing the tube response reduce to ordinary differential equations. 
Interesting and useful engineering Information is obtained from the 
solutions. In addition, It Is possible to see explicitly how the 
assumption of transverse isotropy in a layer affects the differential 
equations governing the behavior of that layer. 

Looking beyond the axlsymmetric problem, the solution to the 
axisymmetric problem can serve as a guide for constructing approximate 
solutions to problems having no closed-form solution. An example of 
such a problem is the determination of the stresses In a layered tube 
with temperature-dependent material properties and subjected to a 
circumferential temperature gradient. Such a condition can exists with 
a tube in space, with the one side heated by the sun and the other side 
exposed to the cold of its own shadow. A significant circumferential 
temperature gradient would exist within the tube. To address such a 
problem, finite-elements, or some other numerical technique can be 
used. Here the use of the principle of complementary virtual work, in 
conjunction with a Ritz approach, is explored. However, here the 


3 


circumferential temperature gradient problem Is not solved. The 
axisymmetric problem Is resolved in an approximate sense, using simple 
functional forms for the stresses. The stresses and deformations 
determined by this approximate approach are compared with the elasticity 
solutions in order to evaluate the use of this approximate approach for 
more complicated problems. 

The paper begins by deriving the governing equations for the 
elasticity solution. Key assumptions and their effect on the governing 
equations are emphasized. The second section presents numerical results 
based on the elasticity solution. The effect of layering arrangement on 
the residual stresses and the axial thermal expansion are illustrated. 
The third section formulates the problem from the point of view of the 
principle of complementary virtual work and a Rltz approximation to the 
stress fields. The fourth section presents a comparison between 
numerical results from the virtual work solution and the elasticity 
solution. The weak and strong points of the approximate solution are 
discussed. Though not directly related to the study of tubes, the fifth 
section presents an interesting sidelight, namely the examination of the 
residual stresses in a solid cylinder. It is demonstrated that the 
theory predicts unbounded stresses at the center of the cylinder if the 
material constants satisfy a certain rather general condition. 



DERIVATION OF THE GOVERNING EQUATIONS 

Figure 1 shows the coordinate system and nomenclature associated 
with the layered tube. The inner radius Is denoted as r^, the outer 
radius r Q , and the radii of the Interfaces r^, r 2 , etc. The axial 
coordinate Is x and the displacement In that direction Is u. The 
circumferential coordinate Is 0 and v denotes circumferential 
displacement. The radial displacement is w and r the radial 
coordinate. The usual notation is used to Identify the components of 
stress and strain, l.e., Is the normal stress in the radial 
direction, y v . is the engineering shear strain in the x-0 plane, etc. 

The temperature of the tube is aT above some reference temperature. 

Here the reference temperature will be the cure temperature of the 
material and aT will be negative. The equations governing a single 
isolated layer will be derived. The solution to the equations will be 
expressed in terms of unknown constants of Integration and the layer 
material properties. There will be one such solution for each layer. 

The constants associated with each layer will be determined by enforcing 
the stress-free boundary conditions on the Inner and outer surface, r^ 
and r Q , and by enforcing certain interface conditions at r^, r 2 , etc. 

In addition, because of the planar nature of the solution, certain 
Integrated conditions will be used in the determination of the 
constants. 

Because the problem is ax i symmetric and each layer Is orthotropic 
In the x-e-r coordinate system, the circumferential component of 
displacement in each layer is zero. Because of axi symmetry, neither of 
the other two components of displacement depend on the circumferential 
coordinate. Furthermore, in the portion of the tubes away from the 
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ends, the radial displacement Is a function only of the radial 
coordinate. With these assumptions, the displacements in each layer 
take the form: 


u = u(x,r) 

v = 0. 

w = w(r) 

The strain-displacement relations then simplify to 

3U , w , dw 

x ax ’ e r * r dr 

v = o • V = A!± • v =o 

Y xe * Y xr ar ' Y er ’ 


( 1 ) 


( 2 ) 


where the total derivative has been used where possible. 

Since each layer is orthotropic, or transversly Isotropic, the 
stress-strain relation for each layer can be expressed in the x-e-r 
coordinate system as 



In this coordinate system thermal expansion effects are purely 
dllatational . The C^'s are elastic constants related to the Young's 
moduli, the Poisson's ratios, and the shear moduli of the material. 

The a's are the coefficients of thermal expansion. Here all properties 
will be assumed to be temperature-independent but the methodology is 
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valid for the temperature-dependent situation. The material constants 
are different for each layer and Appendix 1 gives expressions for 
the In terms of the more familiar engineering constants 
(Ep Eg, .... vp .... ap in the principal material system. 

By eq. 2 and the stress-strain relation. 



e x (r) = Ar + B (6) 

A and B being constants. The second equation requires A to be zero. 
Redefining B to be e°, 

c x (r) = E ° (7) 

and it is seen that the axial strain In a layer is not a function of any 
of the coordinate variables. 

Likewise, with the assumptions of the independence of the stresses 
of x and e, and using eq. 4, the equilibrium equations simplify to 
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( 8 ) 


Tf + r 0 

di i 

xr + _xn = o. 

dr r ’ 


the third one being automatically satisfied. Integrating the second 
equation results In 


xr r 


( 9 ) 


C being a constant. 

The functional form of u(x,r) can be found from the above 
results. Using eq. 7, the first strain-di p placement relation (eq. 2) 
integrates to 


u(x,r) = e°x + f(r) , 


( 10 ) 


f(r) being an arbitrary function of r. From eqs. 2, 3, and 9, 


du df xr _ C 1 


Y xr dr ” dr “ ^ 


55 


C r 
a5 


or 


(ID 


f(r) = f— In r + 0, 
C 55 


D being a constant. With this result 


u(x,r) = e°x + ~ — In r + D. 
C 55 


( 12 ) 


(13) 


By writing the stresses In terms of the displacements, the first 
equilibrium equation leads to a differential equation for w(r), namely 


p r d^W , 1 dw% 7 =r w 1 (r 7? 1 o E 

C 33^ . 2 r dr j ' L 22 2 " r lL 12 " L 13 Jt r ’ 

dr r 


(14) 


where 
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L » (C 13 - C 12 )a x + (C 23 - C 22 )a e + (C 33 - C 23 )a r )6T . (15) 

With the solution to u(x,r) available and by solving for w(r), the 
complete solution Is available for application to either a single-layer 
problem or to a multi pi e-Tayer problem. 


SOLUTION FOR wfr) 

Equation 14 has a solution of the form 


w(r) = A.r* + A 9 r" x + — 

It/ r 


( C 12 “ C 13 ) 


e r e r 


(^33 " 


(16) 


2V 


with 


x 



(17) 


Many times It is assumed that a 1 ayer of fiber-reinforced material Is 
transversely Isotropic In the plane perpendicular to the fibers. In the 
nomenclature of Appendix 1, and In the nomenclature common to fiber- 
reinforced materials, this is the 2-3 plane. While isotropy In that 
plane may not be exactly satisfied, for lack of complete material 
property data, transverse ‘isotropy is often assumed. Referring to 
Appendix 1 for nomenclature, the assumption of transverse isotropy 
Implies 


E 3 = E 2 ; v 13 = v l2* G 23 


2"(t~ J “2 


(18) 


In this problem the assumption Influences the equation governing the 
behavior of a 0° layer. For a transversely Isotropic 0° layer, 


C 33 = C 22 ’ C 13 = C 12* a r = a 0 " 

and the equation for w(r), eq. 14, becomes homogeneous and not dependent 
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on material propei ties. The equation becomes 

p 

d w . L dw w n / 9n \ 

7? + ri‘l {20) 

dr r 

and has solution 

A ? 

w(r) ° A L r + — (21) 

To determine the constants A^, Ag, C, D and e°, boundary. 

Interface, and an Integral condition must be applied. 

BOUNDARY CONDITIONS 

Single Layer 

The conditions for a single layer are that tractions vanish at the 
Inner and outer radii, namely, 

°r< r 1> ‘ °* 'Y< r o> = ° : V (r 1> " 0: V (r o> - 0 

' er < r 1> ” 0i “ 0 (22a ' f) 

and that the net axial force acting on the cross-section of the layer Is 
zero. This axial force condition can be expressed as 

J o rdedr = 0 , (23) 

A x 

where A is the cross-sectional annular area of the layer. The Integral 
can be written more specifically as 

2n J 0 o (r)r dr = 0. (24) 

r 1 

Other resultant equations, such as the net cross-sectional moment and 
the net torsional moment being zero, can be written for the cross- 
section. However, all except eq. 23 are automatically satisfied. As a 
final condition, rigid body displacements of the layer must be 
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suppressed. 

By eq. 4, eqs. 22e and f are automatically satisfied. Application 
of eqs. 22c and d requires C of eqs. 9 and 13 to be zero. Suppressing 
rigid body axial motion forces 0 of eq. 13 to be zero. As a result 

u(x,r) = e°X. (25) 

Applying eqs. 22a and b t and eq. 24 leads to three equations for Al, A2, 
and e°. With these coefficients known, the stresses and displacements 
In a single layer can be determined. 


Multiple Layers 

For multiple layers, the solutions for t (r), u(x,r), and w(r) in 

A I 

eqs. 9, 13 and 16 require a different set of constants for each layer. 
Because the material constants for each layer can be different, the 
value of x will be different for each layer. For the kth layer, 


,00 _ cW 

xr r * 


(26) 


u (k ^(x,r) 


. E o(k) x+ £j> 1nr + D (k) j 

C 55 ; 


(27) 


and 


w ( k )( r ) 


fkl xOO 

a i r + A r r 


(k)„-X (k) + 1^12 


(k) _ *(k) 


c i3 ; ) o< k > 


(C 


33 22 


c^~) Z " 

I 


+ E^r, 


where 


I< k > - <(c< k > - c<k))„(«) t (r (k) . c(k) )o (k) + ( c(k) . c| k 3))4 k )) 4 r, 


and 


00 


Wk) 

u 22 

L 33 


(29) 
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If the layer Is transversely Isotropic, 


v/( k ) (r) 




t (k) 


A}"'r + 


( 30 ) 


'1 • r 

For N layers, there are N e°'s, N C's, N D's, N A1 and N A2's, or 
5N unknown constants. These constants are determined by satisfying the 
traction-free conditions on the inner and outer radii, satisfying 
continuity of tractions at each Interface, satisfying continuity of 
displacements at each Interface, suppressing rigid body motion, and by 
applying an integrated condition on the cross-sectional area. These are 
explained below. 

Since layer 1 Is the inner layer and layer N is the outer layer, 
the traction-free conditions at the Inner and outer radii are 

4‘V,) - 0 ; tjjv,) = 0 ; T<)>(r t ) = 0 (31a-e) 

4 N V 0 > » 0 ! T<|!>(r 0 ) - 0 ; T<|)((r o ) . 0. (31d-f) 

By eq. 4, eqs. 31c and f am automatically satisfied. Equations 31b and 
e lead to, from eq. 26, 

c (0 = C < N > = 0. (32) 
Continuity of the interface tractions can be expressed as 


.00 


< r k> ■ ^ +1) ( r k) ! ' ^. +1) (r k ) 

' 0 r >(r k> * 'er +1)(r k ) • k - ‘-2 N-l. 


(33a-c) 


Again by eq. 4, eq. 33c is automatically satlsfed. Equation 33b, along 
with eq. 32, leads to the conclusion that 


>(k) = 


= 0 


k = 1 ,N. 


(34) 


More will be said of eq. 33a shortly. 
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Continuity of the Interface displacement can be written , 
u (k) (r k ,x) - u< k+1 >(r k ,x) ; v< k >(r k ) ■ v< k+ ’>(r k ) ; 
w (k, ( r k ) - w( k+, )(r k ) , k - 1.* N-i. 


( 35 ) 


Obviously since v = 0 In all layers, continuity of the circumferential 
displacements Is automatically satisfied. Substituting from eq. 27. and 
using eq. 34, eq. 35a results In 

e°< k) x ♦ D< k > . .0(W) X ♦ D < k *D , (36) 

This equation leads to the conclusion that the constant axial strain for 
each layer Is the same for all layers and so the tube-as-a-whole has 
strain e°. This strain is given by 


‘ 0(k) -c° , k = 1,N . 

Also from eq. 36, all D< k ) must be equal and. If axial rigid 
Is eliminated, 


(37) 


body motion 


At this point 


» (k) =0 , k - 1,N 


U^(r,x) = e°> 


and 


(- 8 ) 


(39) 


« (k »(r) - A( k >r k(k) ♦ A^r^’ ♦ ♦ T< k >r 


or 


C (k) p(k) ' 
C 33 C 22 

i(k) 


(40) 


w' ^(r) = A^^r + 1 

' r H 1 r + "p— • (41) 

These expressions for u and w involve N A/s, N A 2 ' S) and e °, or 2N + 1 
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unknowns. The first and fourth of eq. 31, the first of eq. 33, and the 

third of eq. 35 provide 2N equations. The remaining necessary equation 

Is determined by the Integral condition that there Is no net axial force 

acting on the cross section of the tube, l.e., 

J a rdedr, (42) 

A x 

where now A is the annular area of all N layers. Specifically this 
condition is written as 


N r k 

2v z J a^rdr = 0 . (43) 

k=1 Vl 

The complete solution for the displacements, and hence the stresses, for 
the N layers is now available. This whole process can be easily 
automated for an arbitrary number of layers, each with arbitrary but 
orthotropic or transversely isotropic material properties. 

NUMERICAL RESULTS 

To Illustrate some of the thermal effects that occur in layered 
composite tubes, the numerical results for several specific problems 
will be presented. Table 1 lists the material and geometric properties 
used in these numerical examples. 

Effect of Stacking Sequence on Residual Stresses 

Figures 2 and 3 Illustrate the residual axial, circumferential 
(hoop), and radial stresses in two four-layer graphite-epoxy tubes. The 
tube In fig. 2 has a stacking sequence of (0/90/0/90) while the tube of 
fig. 3 has a stacking sequence of (90/0/90/0). Each tube is at a 
temperature 280 °C (500° F) lower than its curing temperature. As far as 
axial stiffness is concerned, these two tubes are practically identical. 
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as are other combinations of two 0° layers and two 90° layers. (Axial 
stiffness Is an important design parameter In many space structures.) 

In the figures the stresses are Illustrated as a function of the 
nondlmenslonal radius p, p = (r - r^) / (r Q - r^) . 

By examining figs. 2 and 3 It can be seen that for both tube 
configurations there are high tensile stresses perpendicular to the 
fiber direction. Also, there are high compressive stresses In the fiber 
direction. The magnitude of the tensile stresses are such that the 
brittle epoxy could easily crack. These high stresses could be expected 
and predicted to some degree of accuracy without resorting to elasticity 
theory. However, It has been found for these particular cross-ply 
tubes, the less refined theories, e.g., classical lamination theory, 
underpredict the compressive stresses In the hoop direction. In 
addition, figs. 2 and 3 show an interesting effect concerning the radial 
stresses that can only be predicted using the more refined theory. 

Figure 2 shows that for the (0/90/0/90) tube the Interface between the 
outer 0° and 90° layers experiences a tensile radial stress. The value 
of the tensile stress is not high compared to the axial and 
circumferential stress components. However, the peak tensile stress 
occurs at an interface between layers, a region that can be weaker in 
tension than the other portions of the tube. Figure 3 shows that by 
interchanging layer order, the interface with the peak stress now 
experiences a compressive stress and the peak tensile stress, which is 
now at another interface, is considerably less. 

Effect of Stacking Sequence on the Coefficient of Thermal Expansion 

As the temperature of the tube is reduced below its cure 
temperature, the thermal expansion effects, coupled with the elastic 
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properties, will produce a change in length of the tube. The 
quantity e° Is a measure of this axial deformation. Recall that e° Is 
the same for all layers. The quantity formed by dividing e° by AT Is 
the axial coefficient of thermal expansion (CTE) of the tube. For 
dimensional stability studies, this expansion coefficient is of great 
interest. As will now be demonstrated, stacking sequence Influences 
this quantity to some degree. 

If all combinations of four-layer tubes with the fibers In two 
layers being axial and the fibers in the two other layers being 
circumferential are considered, six different stacking sequences are 
possible. Table 2 shows the axial CTE's of these six combinations. 
Suprislngly enough, the coefficients vary by 7% among the six tubes. 
While a 7 % variation may seem low, there have been efforts to construct 
tubes with a specific axial CTE. It would seem that stacking arrange- 
ment, as well as fiber orientation, could be used to help achieve this. 

Further Examples 

To provide motivation for the next portion of the paper, and to 
provide other examples of residual effects, a (90/0g/90) tube is 
considered. Such a tube Is a prime candidate for structural 
applications. This tube has a high axial stiffness due to 80% of the 
fibers being In the axial direction. The circumferential layers are 
referred to as skins while the six axial layers are referred to as the 
core. The skins serve to keep the axial layers together. However, 
these circumferential layers contribute significantly to the thermal 
stress state. 

Figure 4 shows the three components of thermally induced stress for 
this skin/core tube. In both the skin and the core the fibers are in 
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compression. Due to the high axial stresses In the skins, the skins 
probably will exhibit circumferential cracks. This does not degrade 
their function of holding the core together. The core Itself could well 
have cracks In It also. The cracks would be radial, running the length 
of the tube. There Is a tensile peak in the radial stress between the 
outer skin and the core. The radial stress between the Inner skin and 
the core Is compressive. 

If the effects of temperature-dependent material properties are 
Important In a particular Instance, the proceeding analysis can be 
used. Temperature dependence of the Q^'s and the o's does not affect 
any of the derivation. However, as mentioned in the INTRODUCTION, the 
effects of circumferential temperature gradients on the stresses could 
be of interest. If the material properties are assumed to be 
independent of temperature, closed-form solutions may be possible, 
depending on the functional form of the gradient. However, If 
temperature-dependent material properties are to be Included with a 
circumferential gradient, closed-form solutions are probably not 
possible, except for some special cases. In general, an approximate 
method would have to be used to determine the stresses. The next 
section explores the use of the principle of complementary virtual work 
in conjunction with a Ritz approximation to the stress fields to study 
thermal stresses in layered tubes. Only the uniform temperature case 
just considered will be studied. However, comparisons of the 
complementary virtual work approach with the elasticity solution will 
serve as a measure of the accuracy of the approximate approach. 

APPROXIMATE METHOD BASED ON COMPLEMENTARY VIRTUAL WORK 

For the problem here, if the volume V is considered as the volume 
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for all N layers, and the surface S Is considered as the Inner surface 
plus the outer surface of the tube, the principle of complementary 
virtual work can be written as 


j(e 60 + e„6a„ + e_6a + y. + y + y 6t )dV 

J v ' x x 0 Q r r 'or or ’xr xr T xe xo' 


+ J (ufiT + vfiT +■ wfiT )dS = 0 
s x 0 r 


m 


For this problem the tractions on the surface, T , T . T , are known, 

x 0 r 

specifically they are zero. Therefore the surface Integral is not 
involved. Implicit In the above equation Is the fact that the stresses 
being varied must satisfy the equilibrium equations and the traction 
boundary and Interface continuity conditions. 

When using the complementary principle It Is more convenient to 
Invert the stress-strain relation of eq. 3 and write (see Appendix 1) 



Approximate solutions to the thermal stress state can be obtained 
from eq. 44 by assuming that the stresses can be expressed in terms of 

I 

specific functional forms and unknown coefficients. Explicitly 
integrating the volume Integral results in Euler equations, In the form 
of linear algebraic equations, which can be solved for the unknown 
coefficients. The nature of the stresses Illustrated In figs. 2-4 
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provides a strong motivation for assuming a simple linear variation of 
stresses within each layer. If circumferential temperature gradients 
were Involved It would make sense to assume the stresses could be 
approximated by the product of function of 0 and linear function of r. 
Thus for an N-layer tube, the stresses are assumed to be of the form 

„(k) . e (k) + c (k) r 

+ b l k)r k = l.N , (46) 

„< k > = a< k > ♦ a( k >r 
r 0 1 

all other stresses being zero. The superscript denotes layer number 

and a£ k ^ , . . . ,cj k ^ are to-be-determined constants. For a N-layer tube 

there are 6N constants. However, to use these In the virtual work 

statement the assumed stresses must satisfy the equations of equilibrium 

and the traction boundary conditions. In addition, for this problem 

there are two further conditions. The first one is that 0 must be 

r 

continuous at the Interfaces between layers. The second one is that the 
axial stresses In the layers must be such that the net axial force is 
zero. 

To satisfy the equilibrium equations, eq. 8, 

a£ k ) = b£ k ) and 2a{ k ^ = bj k >, k = 1,N . (47) 


Since the inner and outer layer (k = 1 and N) are traction free at r * 
r^ and r = r Q , respectively, applying eqs. 31a and d to eq. 46 leads to 


a^ 1 ^ + aPV . = 0 

0 11 

a< N > ♦ aWr . 0 . 
0 1 0 


(48) 


Because o r Is continuous at the Interface between the kth and (k + l)st 
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layers, 


a < k > ♦ a ( k V k . a < M > ♦ a < k+1 >r k . k . l.N-1 


(19) 


The remaining condition on the approximate stresses Is the Integrated 
condition on the cross-section, eq. 42. This leads to 

{ J (c* 1 ) + cj^rjrdr + z J (c^) + c[ k V)rdr 


k=2 r, 


+ / (c< M > + cj^rjrdr} = 0 

r N-l 


(50) 


For a two-layer tube, for example, the stresses that satisfy the above 
conditions are given by 


layer 1 


0 = c< l > + c^r 
°x C 1 C 1 r 


°6 ■ 4 I,(2r • r i> 

° r " a<‘>(r - r,) 


(SI) 


layer 2 


°x - c P ) < r - 8 3> - a l c l l> - 6 2 c , 


( 1 ) 


°e ■ ‘i 0 ( Fp-f 0 >' 2r - r o> < 52 > 

The quantities s^, f^, and e 3 are known constants involving the tube 
geometry. They are defined in Appendix 2. There are four constants, 
aj 1 ), c^, cl 1 ), cp) , that determine the stress state. For a three 
layer tube there are seven constants. 
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Substituting the stresses Into the strains and the strains Into eq. 
44, and also substituting the first variations of the stresses Into eq. 
44, Integration of eq. 44 leads to an equation of the form 



The equations 



(53) 


T i a 0 , 1=1,4 

lead to solutions for aj 1 *, c^, cj 1 ) and c{ 2 ). 


(54) 


NUMERICAL EXAMPLES USING APPROXIMATE METHOD 
Figures 5-7 show examples of the predictions of the approximate 
method. The exact elasticity solution Is Included In each example - as 
a comparison. The examples depict situations designed to Illustrate the 
level of accuracy of the approximate method rather than to Illustrate 
results of engineering significance. 

Figure 5 shows the three components of stress for a (0 3 /90) tube. 
Though this really is a four-layer tube. It Is modeled as a two-layer 
tube, the 0° layer being three times thicker than the 90° layer. 

Despite this lumping of the layers, the comparison between the 
approximate approach and the exact solution Is quite good. The 
correlation seen here Is about as good as has been observed In these 
studies. As has been seen with the elasticity solution, the peak radial 
stresses generally occur at the Interface between layers of dissimilar 
orientation. It has been found that the approximate method will accurately 
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predict the peak value of radial stress at the Interfaces even though at other 
locations the predictions for this same stress are not as good. 

Figure 6 shows the results for a ( 9 O 3 /O) tube. Again the three 90° 
layers are lumped together and are represented as one layer. While the 
axial stress and the peak radial Interface stress are accurately 
predicted, the circumferential stresses In the 90° layers are not 
accurately predicted. Similarly, the predictions of the radial stress 
In the 90° layers Is poor. The correlation seen here Is about as poor 
as observed In this study. It has been found In these studies that the 
thicker the grouping of layers with circumferential fibers, the less 
accurate Is the approximate method. The principle of complementary 
virtual work with linear stress approximates tends to underpredict the 
peak circumferential stress. For these residual stress states, with a 
negative aT, the inaccurate circumferential stresses represent 
compression in the fiber direction. Since compression in the fiber 
direction is not a concern In these tubes, this error, at least In this 
study, is tolerable. Of course, when using the approximate method, a 
further discretization of the grouped layers helps Improve the 
prediction. Figure 7 shows the effects of splitting the one three-layer 
grouping In the (90 3 /0) In half and representing It with two layers. 

The tube Is then represented by three layers. The principle of 
complementary virtual work uses the extra degrees of freedom quite 
efficiently. With just this refinement It Is easy to predict the 
results of addition discretization of the group of 90° 's. 

Finally, fig. B shows the comparison between the approximate method 
and the theory of elasticity for a tube of engineering interest, the 
(90/0g/90) tube examlneu before. It is clear that the fairly crude 


22 


approximation to the stress fields provides a good estimate of the 
stress state. It Is expected that when a circumferential temperature 
gradient and temperature-dependent material properties are considered 
for this type of tube, the linear approximation to the stress fields 
should be sufficiently accurate for parameter studies. 

DEFORMAT I ONS PREDICTED USING THE PRINCIPLE Of COMPLEMENTARY VIRTUAL WORK 
While the principle of complementary virtual work Is Ideal If 
concern Is solely with the stress state, It Is of Interest to study the 
deformations predicted by the principle and compare the results with the 
predictions of elasticity. 

By using the assumed stress fields In the stress-strain relation of 
eq. 45, the strains for the kth layer can be written In the form 

E ( k ) = + a^aT + B^r 

X X 

e (k) a + aj^aT + D^r (55) 

0 0 

e ( k ) = + a^AT + F^ k V 

r r 

There are no shear strains. In the above the quantities A^ , . . . ,F^ k ^ 
are functions of the layer compliance, S^, and the constants In the 
assumed form of the stresses, a^, ..., c[ k ^, for that layer. Not 
surprisingly, the above form of the strains do not satisfy the 
compatibility equations. Therefore any predictions regarding 
deformations must be viewed with caution. However, it Is Instructive to 
examine the predicted axial CTE for several of the cases shown in Table 
1. Table 3 Indicates the values of the axial CTE's for four of the six 
tubes of Table 1. Only these four were examined because the approximate 
analysis was limited to tubes with a maximum of three layers. The 
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term a x In Table 3 is the free thermal expansion of each particular 
layer In the axial direction. They are taken from Table 1. The term 
labeled 'approximate axial CTE ' Is the quantity { (A/aT) + a x ) and from 
eq. 55 Is the portion of e that does not vary with r. The terms 

A 

A/aT and B/aT represent the effect of the other layers on the free 
axial expansion of the particular layer. It Is clear that the B/aT term 
contributes very little to e This might be expected since the 

A 

elasticity solution predicts e to be Independent of r. The fourth 
column of Table 3 Is the axial CTE predicted by the theory of 
elasticity. (Recall, the theory of elasticity predicts that the axial 
CTE Is the same in all layers.) The agreement between theory and the 
approximate method Is quite good. With the breakdown In Table 3 It Is 
Interesting to see how A/aT contributes to the overall axial expansion 
of a particular layer. For example. If a Is large, A/aT Is negative. 

A 

This means that the other layers must restrain the large free thermal 
expansion value so the deformation of the layer Is, to some degree, 
compatible with the deformations of the other layers. Since 
compatibility Is not satisfied, and since the approximate solution makes 
no explicit statements regarding continuity of Interface displacements, 
the results of Table 3 are remarkable. 

THERMAL STRESSES IN A SOLID CYLINDER 
As a final note s consider thr case of a solid and unlayered 
orthotropic cylinder. Equations 9, 13, and 16 govern the behavior of 
the cylinder. In this case there is no inner radius at which to apply 
the boundary conditions. Rather, statements must be made reyarding 
boundedness at the origin, as well as satisfying the boundary conditions 
on the outer surface, and the integral condition. Specifically, 


24 




( 56 ) 


o (r ) = 0 : t (r ) = 0, 
o' * xr v o' * 

and w(r) and u (x,r) remain bounded as r + 0, 

l.e., w(r), u(x,r) < », r «■ 0. (57) 

Also, 

J o rdrde. (58) 

A x 

Since x > 0, the boundedness condition on w(r) of eq. 16 requires 

A 2 « 0 . (59) 

The condition on the shear stress t at r = r and the suppression of 

Al ^ 

axia' 1 rigid body translation requires 

C = 0 (eq. 9) , 

D = 0 (eq. 13). (60) 

The constants Aj and e° can be determined from the condition of the 
vanishing of the normal stress at r = r Q and the Integral condition. In 
general, neither A^ nor e° are zero. 

For the solid cylinder the nonzero stresses are proportional to 
r x ~*. If C 22 > Cg^, then x Is greater than unity. At r * 0 all three 
components of stress are zero. On the other hand, If < C 33 , then 
x is less than unity and the stresses are unbounded at r = 0. The 
case C 22 < corresponds to the situation when the fibers In the 
cylinder are oriented radially. It may be difficult to manufacture such 
a cylinder and have the material properties remain constant with 
radius. Nonetheless, for this case the tube would be much stiffer in 
the radial direction than in the circumferential direction. Figure 9 
shows the three components of stress for the case of using the material 
properties of Table 1 and assuming the fibers are oriented radially. 

The value of x is 0.298. It should be mentioned that the value of x is 
a material property and is independent of the fact that the problem is a 
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near the 


thermal stress problem. The stresses are proportional to r*"^ 
origin for the case of mechanical loads, e.g., an axial load, being 
applied to the cylinder. 

While the Infinite stress Is physically Impossible, the coefficient 
of the Infinite stress Is of Interest. The coefficient is an indication 
of the severity of the stress gradient near the centerline of the 
cylinder. Unlike the numerical value of \, this coefficient Is a 
function of whether or not the cylinder is loaded thermally or 
mechanically. The coefficient also depends on the radius of the 
cylinder. In fracture mechanics this coefficient is related to the 
stress intensity factor. Table 4 presents the stress Intensity factors 
for the three nonzero components of stress. It is clear from the table 
that the radial stress Is the most intense as the centerline of the tube 
is approached. Since the three intensity factors are positive, cooling 
the tube results in high compressive stresses near the centerline. This 
situation is more favorable than having high tensile stresses, 
particularly in the matrix dlre^Jon. 

CONCLUSIONS 

From the results presented, it can be concluded that for cross-ply 
composite tubes cured at an elevated temperature, residual thermal 
stresses can be of significant nagnitude. The peak stresses, as well as 
the axial CTE, are influenced by the stacking arrangement. It appears 
that the principle of complementary virtual work, in conjunction with a 
Ritz approach, can be used to estimate the thermally-induced stresses. 
The approach also can be used to predict the axial deformations with a 
remarkable degree of accuracy. This, despite the lack of the 
approximate solutions satisfying the compatibility equations. When 
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considering approximate solutions to thermal stress problems It would 
seem that stress-based approaches, such as the principle of 
complementary virtual work, are preferrable over displacement-based 
approaches. For thermal stress problems, displacements need not be 
considered and so there Is no need to approximate variables which must 
be differentiated to determine the real quantity of Interest. Finally, 
in the idealization of linear elasticity, It appears solid orthotropic 
cylinders can experience high stress gradlants at their centers. 
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TABLE 1 


Material and Geometric Properties 

E 1 = 146.8 GPa; Eg - 9.929 GPa; E 3 = 9.101 GPa psl 

v 12 = v 13 ~ v 23 ~ 

0l = -0.0774 x 10" 6 /°C; * 2 - 33.66 x 10" 6 /°C = a 3 

r.j = 6.35 mm, layer thickness = 0.127mm 
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TABLE 2 


Axial Coefficients of Thermal Expansion 
For Cross-Ply Tubes 

Tube axial CTE (per °C) 


(90 2 /0 2 ) 

2.468 x 10 -6 

(90/0/90/0) 

2.538 x 10" 6 

(0/90 2 /0) 

2.601 x 10“ 6 

(90/0 2 /90) 

2.612 x 10" 6 

(0/90/0/90) 

2.675 x 10" 6 

(0 2 /90 2 ) 

2.740 x 10" 6 
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TABLE 4 


Stress Intensity Factors For Solid Cylinder 
$22 < C 33 ^ 


Stress Component 


Coefficient of r x ~* 


0.0848 

0.1303 

0.4366 


* units are such that stress is in MPa/°C increase in temperature 
above the reference temperature. 
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Fig. 1 Geometry and nomenclature of the tube cross-section. 
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Fig. 5 


elasticity 


Comparison between approximate 
and exact solutions for a 
(0 3 /90) tube. 






axial stress. MPa 










hoop stress, MPa axial stress, MPa 


100 




radial stress, MPa 



approximate 


elasticity 


Fig. 8 Comparison between approximate 
and exact solutions for a 
(90/0 6 /90) tube. 
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Appendix 1 

Details of the Stress-Strain Relation 


Figure A1 shows a segment of a cylindrical lamina and the principal 
material 1-2-3 and the cylindrical x-e-r coordinate systems. In the 
1-2-3 system: 
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i 

if 

The Inverse of this relation is: » 
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and eq. A-3 transforms Into 



In the above 


S 1 1 - m + m n ( 2S + Sgg) + ^22 

^ i 2 = $2\ = ( m4 + n4 ^ S 12 + + S 22 " S 66 ^ 

^13 = ^31 = S 23 + m S 13' 1 

S 16 ‘ ron (m 2 (2S^ - 2$ 12 - s gg) + n ( 2S i2 “ 2S 22 + S 66^ 
^22 = n ^11 + m n ( 2 ^12 + ^ 66 ^ + m ^22 

^23 = ^32 = S 23 + n S 13^ 

^26 = mn ( 2 ^12 “ ^ 22 + ^ 66 ^ + n ^11 “ 2 ^12 " ^ 66 ^ ^ 
^33 = S 33 

^36 = $63 = 2mn ( s l3 - S 23^ 

- 2 2 
S 44 = m S 44 + n s 55 

^45 = ^54 = mn ( S 55 _ W 

^55 = n S 44 + m S 55 

Sgg = 4m n (S^ “ ^^12 + **22^ + ^66 ^ n ) 
m = cos(|) n = s i na> . 
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Also 

a x = m ai + n a 2 ; a 0 = n 2 ai + l/c^, a r = a 3 
°or “ 0 ‘ a xr’ a xe = 2mn ( a i - a 2^ ♦ 


C n = m Z, C u + 2m 2 n 2 (C 12 + 2C gg ) + n 4 C 22 
C 12 = n 2 m 2 (C u + C 22 - 4C gg ) + ( n 4 + m 4 )C 12 
C 13 = m Cjj + n C 23 

C 16 - mn[m 2 (C n - C 12 - 2C gg ) + n 2 (C 12 - C 22 + 2C gg ) j 
C 22 = n 4 C u + 2nV(C 12 + 2C gg ) + m 4 C 22 
C 23 - n 2 C 13 + m 2 C 23 

C 26 - ran[n 2 (C 11 - Z n - 2C gg ) + m 2 (C 12 - C 22 + 2C gg ) ] 
C 33 = C 33 

C 3g - mn(C 13 - C 23 ) 

^44 = m ^ C 44 + n ^ C 55 
C 4g = mn(C 5g - C 44 ) 

^55 = n ^ C 44 + m ^ C 55 

^66 = "^ll * 2C 12 + C 22 ) + C gg (m 2 - n 2 ) 2 . 
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Appendix 2 

Definition of Constants In Approximate Method 


5 1 = V H 2 ; 0 2 = H l /H 2 ; 0 3 = G 2 /H 2 


F 2 = ( r i - rJ)/3; Hj = (rj - rf)/2; Gp = (r^ - 


(r? - r?W^> u - /-2 


rf)/3 


H ? * < r o - rf)/2 • 


= inner tube radius 
r 0 = inner tube radius 
r l - radius to interface between layers. 
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